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Abstract. WSRT observations have been used to investigate the presence of neutral hydrogen in extremely young 
radio galaxies. These objects were selected from a sample of High- Frequency Peakers (HFPs). We detect 2 of 
the 6 observed galaxies confirming previous detection of H I in these objects. In the case of OQ 208 - for which 
discrepant results were available - we confirm the presence of a broad (~ 1800 km s _1 ), blue-shifted and shallow H I 
absorption. No significant changes in the Hi profile have been found between the two epochs of the observations. 
The intriguing result is that the derived H I column densities and upper limits obtained for the most compact 
sources, do not follow the inverse correlation between the column density and the linear size found for CSS/GPS 
sources. This can be explained - assuming the gas is already in a torus/disk structure - by a combination of the 
orientation and the extreme compactness of the sources. In this scenario, our line of sight to the source would 
intersect the torus in its inner region with low optical depth due to high spin and kinetic temperatures. There 
is no evidence, with the exception of OQ 208, of unsettled, high column density gas still enshrouding the young 
radio sources. This can be due to the low filling factor of such a medium. 
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1. Introduction 

A substantial amount of gas and dust is often found in the 
central regions of galaxies harbouring an Active Galactic 
Nucleus (AGN). This gas provides important information 
on the physical conditions and processes in these near- 
nuclear regions (Morganti et al. 2004a; Peck et al. 1999). A 
component of this gas, the atomic hydrogen, can be stud- 
ied in radio-loud objects via absorption detected against 
the strong continuum source. Since years, a large number 
of detailed studies has been done (see e.g. Heckman et 
al. 1983; van Gorkom et al. 1989; Morganti et al. 2001, 
2004a, 2005; Vermeulen et al. 2003; Pihlstrom et al. 2003) 
and they have shown that the neutral hydrogen can be 
associated with different structures, from circumnuclear 
tori to fast outflows. These results have been crucial for 
constructing our present picture of the nuclear structure 
of radio loud AGN. 

The medium enshrouding the nuclear regions also 
plays an important role in the growth and evolution of the 
radio source (Labiano et al. 2006). The presence of signif- 
icant amount of gas in young radio galaxies is supported 
by a larger incidence of Hi absorption in these objects 
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(Vermeulen et al. 2003; Pihlstrom et al. 2003) compared 
to what typically found in old and larger radio sources 
(van Gorkom et al. 1989; Morganti et al. 2001). This dense 
medium is likely the result of the merger that triggered the 
radio source (Morganti et al. 2004b). Thus, in the initial 
phase of the radio source, the radio jet likely has to make 
its way through a relatively dense medium that surrounds 
the active nucleus. 

The radio galaxies classified as "Compact-Steep 
Spectrum" and "GHz-Peaked Spectrum" are believed to 
represent the early stages of the radio activity (Fanti et al. 
1995; Readhead et al. 1996; Snellen et al. 2000). These ra- 
dio sources are intrinsically small (sub-galactic size, i.e. < 
15 kpc) and bright (Pi.4GHz > 10 25 W/Hz), characterised 
by a convex radio spectrum, peaking at frequencies be- 
tween 100 MHz and a few GHz. Their radio morphologies 
appear to be the scaled-down version of powerful edge- 
brightened radio galaxies, with luminous mini-lobes (~ 
0.1 up to few kpc) and weaker jets and cores. Both kine- 
matic (Polatidis & Conway 2003), and spectral (Murgia 
2003) studies strongly support the youth scenario, indi- 
cating ages of 10 3 - 10 5 years. 

Given the relatively high detection rate of H I absorp- 
tion in these objects (Vermeulen et al. 2003; Pihlstrom 
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et al. 2003), it is of particular interest to investigate the 
characteristics of the nuclear interstellar medium (ISM) in 
even younger radio sources. In the youth scenario, the anti- 
correlation between the turnover frequency and the linear 
size (O'Dea 1998), which is indicative of the age, suggests 
that the youngest sources have the highest turnover fre- 
quency. Therefore, the "High-Frequency Peaker" (HFP) 
radio sources, characterised by the same properties of 
CSS/GPS, but with the spectral turnover occurring at 
frequencies higher than 5 GHz, are good candidates to be 
newly born radio sources, with ages of about 10 2 -10 3 years 
(Dallacasa 2003). 

This paper reports on the result of observations search- 
ing for H I absorption in a sample of 6 HFP radio sources, 
selected from the Dallacasa et al. (2000) bright HFP 
sample, and suitable to be observed at the Westerbork 
Synthesis Radio Telescope (WSRT). 

2. The sample 

The bright HFP sample (Dallacasa et al. 2000) has been 
constructed by cross-correlating the Green Bank Survey 
(87GB) at 4.9 GHz and the NRAO VLA Sky Survey 
(NVSS) at 1.4 GHz. Only the sources with a flux den- 
sity brighter than 300 mJy at 4.9 GHz and with a ris- 
ing spectrum steeper than a = —0.5 (Soc v~ a ; i.e. to 
avoid turnover frequencies < 5 GHz) have been selected 
as HFP candidates. To tackle the contamination of vari- 
able sources, like blazars, all the candidates have been 
simultaneously observed with the VLA at eight different 
frequencies from 1.4 to 22 GHz. The final sample consists 
of 55 sources (~ 3% of the sources exceeding 300 mJy in 
the 87GB): 11 galaxies with redshifts ranging from 0.02 to 
0.67; 33 quasars with higher redshifts, typically between 
0.9 and 3.5; 5 BL Lac objects, while 6 sources still lack 
optical identification (Dallacasa et al. 2002; Dallacasa et 
al., in preparation). Since the selection of these sources 
is based on their simultaneous radio spectra at a single 
epoch, there is still some contamination by beamed ra- 
dio sources, whose emission is temporarily dominated by 
a self-absorbed jet component. Further multi-frequency 
VLA observations to exclude the long-term variability 
(Tinti et al. 2005) and pc-scale resolution morphological 
analysis (Orienti et al. 2006), have been carried out to 
confirm genuine and young HFPs. 

As shown by Orienti et al. (2006), there is a clear 
segregation in radio morphology between galaxies and 
quasars. The majority of galaxies (9/11) show a Compact- 
Symmetric-Object(CSO)-like morphology, while 28/33 
quasars display "Core- Jet" or unresolved structures. This 
is consistent with the idea that the HFP spectrum in 
galaxies and quasars originates in different regions: mini- 
lobes and/or hot-spots in galaxies, compact regions re- 
lated to the core and the jet base in quasars, in agreement 
with what found in GPS (Stanghcllini et al. 2005) and 
bright CSS samples (Dallacasa et al. 1995; Fanti et al. 
1990). 

Since the radio sources hosted in galaxies likely represent 



a population of young radio sources, six out of the eleven 
HFP galaxies, have been observed to search for Hi ab- 
sorption. The selected objects have recessional velocities 
such that the redshiftcd Hi line is shifted to frequencies 
that are known to be relatively free of radio interference 
(RFI) at the WSRT. Their redshift ranges from 0.02 and 
0.67, and the flux densities at 1.4 GHz exceed 90 mJy. 

For one of the five remaining HFP galaxies, WSRT 
observations were precluded, since a well known RFI 
(GSM) affects the frequency band, while for the other four 
sources, the optical redshift is not available. 

3. WSRT observations and data reduction 

The WSRT observations have been carried out in differ- 
ent runs from 29 January to 12 September 2005 in dual 
orthogonal polarisation mode. Two sources in the sample 
were already known to have H I absorption: from Carilli et 
al. (1998) in the case of J01 11+3906 and Vermeulen et al. 
(2003) and Morganti et al. (2005) in the case of OQ 208. 
In these two sources, the goal of the new observations was 
to investigate the possible presence of variability in the H I 
absorption and/or a broad component in their absorption 
spectrum. In particular, in the case of OQ 208, the two 
previous observations report different characteristics (i.e. 
width) of the absorption. 

Of the six observed galaxies, those with z < 0.3 have 
been observed for 6 hours each, using the L-band receiver, 
with 1024 spectral channels covering a 20 MHz wide ob- 
serving band. The sources with redshift higher than 0.4 
have been observed for 12 hours each, using the UHF-high 
band, with 1024 spectral channels. The 20 MHz band, al- 
lowed to cover a velocity range around the velocity cen- 
troid of approximately ± 2200 km s -1 at z = 0.07, and 
approximately ± 3500 km s _1 at z=0.6. Only the source 
J0003+2129 has been observed with a 10 MHz wide band 
(<~ ± 1500 km s _1 ), since known RFI would have affected 
a wider band. Table 1 summarizes the observational and 
physical parameters of the sources. 

The data reduction was carried out using the MIRIAD 
package. The quasars 3C 286, 3C 48 and 3C 147 were ob- 
served to calibrate the absolute flux density scale and 
the bandpass. Data were inspected for time-limited and 
baseline specific interference, and bad data were removed 
before solving for the gain- and band-pass calibrations. 
The two orthogonal linear polarisations were added to- 
gether to improve the signal-to-noise ratio of the final 
spectrum. The continuum subtraction was done by us- 
ing a linear fit through the line-free channels (i.e. all the 
channels outside a range of ± 1000 km s _1 from the cen- 
troid frequency) of each visibility record and subtracting 
this fit from all the frequency channels (using the MIRIAD 
task "UVLIN"). Particular attention has been given to 
the case of OQ 208 where shallow absorption has been 
detected (see below). After several trials, the final con- 
tinuum subtraction was obtained using as line-free those 
channels outside the range —1700/ + 600 from the centroid 
frequency. The baselines of the spectra obtained after the 
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Table 1. Physical and observational parameters of the 6 HFP galaxies observed with WSRT. Columns 1, 2: source names; 
Column 3: projected linear sizes (Orienti et al. 2006); Column 4: optical redshift; Column 5: central frequency; Column 6: 
channel resolution; Column 7: la noise level in the line cube; Column 8: continuum flux density taken from our WSRT data at 
the observed frequency; Column 9: peak flux density of the absorption line, measured on the spectral image; Column 10: optical 
depth; Column 11: the width of the Hi absorption line: for J0111+3906 the FWHM is given, in the case of OQ208 we give the 
FWZI, due to the complexity of the line profile. Column 12: Hi column density derived from Nm = 1.82 x 10 18 T S pinT peak AV A 
cm -2 , a T S pi n of 100 K has been assumed; Column 13: the redshift of the peak Hi absorption. The line flux density, the optical 
depth and the Hi column density upper limits have been computed assuming the 2a noise level, a line width of 100 km s _1 and 
T sp in of 100 K, as in Vermeulen et al. (2003). a: For the source J0111+3906, the projected linear size is taken from Owsianik 
et al. (1998). 



continuum subtraction result flat, indicating that no broad 
and shallow line features have been included. 

After the continuum subtraction, line cubes were pro- 
duced using uniform weighting, Hanning smoothed and 
inspected. For the two sources where H I absorption was 
detected (J0111+3906 and OQ 208), the spectra obtained 
at the location of the peak of the radio continuum are 
shown in Figs. 1 and 2. For every source, the r.m.s. noise 
was estimated from the line cube and the values are given 
in Table 1. 

In order to make a reliable comparison between our 
values and what found in CSS/GPS samples, the upper 
limits to the optical depth and the Hi column density 
have been computed following Vermeulen et al. (2003), 
i.e. assuming the 2cr noise level, a line width of 100 km 
s" 1 and T spin of 100 K (sec Tab. 1). 

Finally, we have used the line-free channels to pro- 
duce a continuum image of each target. The continuum 
emission is always unresolved at the resolution of the 
WSRT (typically about 20"), with the only exception of 
J01 11+3906, known to possess a relatively weak compo- 
nent on such angular scale (Baum et al. 1990). The con- 
tinuum emission derived from our data is slightly resolved 
to the East, consistent with the large-scale emission seen 
by Baum et al. (1990). 

4. Results 

Of the six galaxies observed, we detect H I absorption only 
in the two sources (J0111+3906 and OQ208) in which 
the presence of absorption was already reported by previ- 
ous studies. In addition to this, we find H I in emission in 



two companion galaxies of J0655+4100. We discuss below 
these three objects in more detail. 

4.1. The source J01 11+3906 

The radio source J0111+3906 is associated with a narrow 
emission line galaxy at z — 0.668 (Carilli et al. 1998). 
The source is classified as a Compact Symmetric Object 
(CSO), due to its VLBI-scale morphology (Taylor et al. 
1996), but it also shows an extended emission on the kpc- 
scale, interpreted in terms of recurrency of the radio ac- 
tivity (Baum et al. 1990). 

Hi absorption was already reported by Carilli et al. 
(1998), based on narrow-band WSRT observations. We 
detect a strong Hi absorption line peaking at ^m.pcak = 
851.03 MHz (z Hl , P cak = 0.6687, see Fig. 1). The Hi ab- 
sorption is only detected against the peak of the contin- 
uum, coincident, therefore, with the VLBI-scale structure. 
The line width is FWHM = 100 km s _1 and the optical 
depth Tp Ca k = 0.437. With these values, assuming a spin 
temperature of 100 K, we obtain an Hi column density 
of 8.05 x 10 21 cm~ 2 . These values are in good agreement 
with those found by Carilli et al. (1998). 
No broad absorption (i.e. with FWHM > 100 km s _1 ), 
with optical depth higher than Ti a ~ 0.03, has been de- 
tected. 

4.2. The source OQ208 

In the radio source OQ 208 we detect a shallow (T pea k 
= 0.005) and very broad (FWZI - 1800 km s" 1 ) Hi 
absorption (Fig. 2a). This absorption extends from <~ 
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Fig. 1. The H I absorption profile detected in J0111+3906. The 
velocity is in the optical heliocentric convention. 



21650 to ~ 23450 km s _1 , therefore mostly blue-shifted 
compared to the systemic velocity 22957 km s _1 from 
Gallimore et al. (1999). The line profile has a peak located 
at ^ H i, peak = 1318.6 MHz (z H i, P cak = 0.0773). Assuming 
a T sp i n of 100 K, we obtain an Hi column density of 8.0 
x 10 20 cm" 2 . 

The Hi absorption in OQ 208 was already detected by 
both Vcrmeulen et al. (2003) and Morganti et al. (2005) 
although with different results. Vermeulen et al. (2003) 
found an H I absorption centered on the same frequency 
(y = 1318.9 MHz) as in our observations but with a width 
of 256 km s~ 1 only. The likely reason for the discrepancy 
with our results is the fact that the maximum available 
bandwidth at the WSRT at the time of the survey of 
Vermeulen et al. (2003) was not large enough to allow the 
detection of such broad features. More recently, the pres- 
ence of a broad Hi absorption in OQ 208 was reported 
by Morganti et al. (2005) and interpreted as a fast out- 
flow of neutral hydrogen. Fig. 2b shows the comparison 
between the H I profile derived from the observations pre- 
sented here and what was found by Morganti et al. (2005). 
In Fig. 2c, the residuals obtained after subtracting the pro- 
files taken in two different epochs are shown. The figure 
shows that there is no significant evidence of changes in 
the H I profile corresponding to an optical depth upper 
limit of t ~ 0.002. 

It is worth mentioning that the radio source OQ 208 is 
associated with a broad emission line galaxy at z=0. 07658. 
The radio structure is characterised by two compact lobes, 
separated by about 10 pc, highly asymmetric in terms of 
flux density ratio (Stanghcllini et al. 1997). In this source, 
the shape of the radio spectrum could either arise from 
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Fig. 2. a: The Hi absorption profile detected in OQ208 from 
our new WSRT observations. The velocity is in the optical he- 
liocentric convention. Strong RFI is present in the line profile. 
6: The two-epoch line profiles of OQ 208 one superimposed to 
the other; the continuous line represents data from Morganti 
et al. (2005), while the dotted line represents our new data, 
c: the residuals obtained after subtracting the two-epoch line 
profiles. 



synchrotron self-absorption or free-free absorption by an 
external absorber embedding the lobes (Kameno et al. 
2000). Furthermore, X-ray observations (Guainazzi et al. 
2004) have led to the discovery of a Compton-thick ob- 
scured AGN where the column density of the absorber is 
10 21 cm~ 2 . As suggested by Guainazzi et al. (2004), in 
this source the jets are likely piercing their way through 
this Compton-thick medium which is pervading the nu- 
clear environment. The outflow detected in Hi would be 
another indication of this process. 
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Fig. 3. Hi total intensity contours superimposed onto a DSS 
optical image centered on the radio galaxy J0655+4100. Two 
nearby galaxies are clearly detected. Contours levels are: 2.35, 
3.50, 5.68 and 7.51 x 10 20 atoms cm" 2 . 



4.3. The source J0655+4100 

Although no evidence of neutral hydrogen gas has been 
found in the radio galaxy J0655+4100, an elliptical galaxy 
at z = 0.02156 (V sys = 6464 km s _1 ; Marcha et al. 1996), 
we find instead, H I in emission in two nearby galaxies. 
One companion is the spiral galaxy UGC03593 where, 
centered at 6694 km s^ 1 , ~2x 10 9 M Q of Hi have been 
detected. This galaxy is located 65 kpc South-East of 
J0655+4100. The second companion displays a systemic 
velocity of 6334 km s _1 and is located to the North- West 
at 55 kpc from the HFP source. In the Digitized Sky 
Survey (DSS) image it is associated with a faint galaxy for 
which we find an Hi mass of ~ 2x 10 8 Mq. This galaxy 
has not been previously catalogued 1 . The two galaxies 
form a physical compact group with J0655+4100 (Fig. 
3) supporting the idea that young radio sources reside 
in groups, as suggested by other works (Stanghellini et 
al. 1993), based on the excess of galaxy density in the 
optical images of GPS radio sources. 



5. Discussion 

The detection in H I absorption of only two of the six ob- 
served galaxies is somewhat surprising. Analysing a sam- 
ple of 41 CSS/GPS radio sources, Pihlstrom et al. (2003) 
have found an inverse correlation between the source lin- 
ear size and the H I column density: smaller sources (< 0.5 
kpc) have larger H I column density than the larger sources 
(> 0.5 kpc). Following this correlation, one would have ex- 
pected very high Nhi in the smallest sources, like HFPs. 



1 Based on results from the NASA/IPAC Extragalactic 
Database (NED). 



In Fig. 4 we show the inverse correlation found by 
Pihlstrom et al. (2003) with the addition of the values 
obtained from the sources presented in this paper. The 
sensitivity reached by the new observations, enables us 
to set tight upper limits to the Nhi for the sources with 
no Hi detection. Interestingly, in our sample, the smallest 
sources (with size ~ 10 pc) appear to deviate from the ex- 
trapolation of the trend found by Pihlstrom et al. (2003). 
Given the limited number of objects, we cannot establish 
whether for very compact radio sources a break down of 
the correlation occurs. If we consider that 15 out of 41 
sources of Pihlstrom et al. (2003) sample have H I column 
density below 2xl0 20 cm -2 , the probability to find two 
sources randomly chosen both below such limit is ~ 13%. 
This means that the small column density is possibly a 
feature typical of extremely young radio galaxies. We also 
point out that the correlation from Pihlstrom et al. (2003) 
has a large scatter, and upper limits to the column density 
are present at all linear sizes. With all this in mind, it is 
interesting to consider the possible origin of the lack of H I 
absorption found in J0655+4100 and J1623+5524. 

The trend observed by Pihlstrom et al. (2003) has been 
explained with both a spherical and axi-symmctric gas dis- 
tribution, with a radial power law density profile, as well 
as a disk model. Indeed, a torus- like distribution of neutral 
gas has been found to be consistent with the observations 
in a number of GPS sources (i.e. 4C 31.04 by Conwayl996; 
1946+708 by Peck et al. 1999). In this scenario, the ab- 
sence of high H I column density in very compact sources 
can be explained by both the orientation and the extreme 
compactness of the sources. Within the framework of the 
AGN unification scheme, the central engine is surrounded 
by a disk of ionised gas and shielded by an obscuring torus 
of atomic and molecular gas (i.e. Urry & Padovani 1995). 
Evidence of ionised disks with radii between a few tens and 
one hundred pc has been claimed in this class of sources 
(Kameno et al. 2000), as responsible of free-free absorp- 
tion in the optically thick region of the radio spectrum. 

The fundamental ingredient to trace H I absorption is 
the presence of background flux density. Given the ex- 
tremely small linear sizes of the HFP sources (see Fig. 4 
and Tab. 1), the line of sight likely passes mainly through 
the ionised region, piercing the torus only in its very inner 
part, if we are looking at the receding lobe. This could be 
the case if the torus axis is aligned to the radio jet, and its 
half-opening angle is, for example, 8 ~ 20° (as illustrated 
in Fig. 5), consistent with the constraint 9 < 30° — 40° 
by Granato et al. (1997), and an orientation of the line of 
sight between 45° and 70° (with respect to the radio jet). 

It should be noted that the gas closer to central en- 
gine is hotter, both in terms of kinetic and spin temper- 
ature, therefore a longer path length would not necessary 
imply a higher optical depth. Moreover, recent work by 
Hatziminaoglou et al. (2006) focused on the properties of 
the dusty tori in AGN, has shown the possibility of low 
optical depth tori. Therefore, in this scenario, it is possi- 
ble that for these sources of a few pc in size (or even less), 
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Fig. 4. Absorbed Hi column density versus projected linear 
sizes. The HFP galaxies do not seem to follow the trend found 
for CSS/GPS sources. The CSS/GPS values have been taken 
from Pihlstrom et al. (2003). 



we are missing the Hi absorption due to both sensitivity 
limitations and orientation effects. 

As described in Sect. 1, in addition to the presence of 
circumnuclear tori, in such young radio sources one could 
have expected to detect also the presence of unsettled gas, 
i.e. gas not necessarily located in a circumnuclear torus or 
disk-like structure but still enshrouding the radio source. 
We do not find evidence for this. Again, this could be 
partly due to the very small emitting area of the back- 
ground radio source, that decreases the probability to de- 
tect the absorption. Clouds indicating the presence of a 
rich and clumpy interstellar medium in the centre of the 
radio source have been found as in the case of 4C 12.50 
(Morganti et al. 2004b), 3C49 and 3C 286.3 (Labiano et 
al. 2006). In 4C 12.50 a cloud of - 20 x 60 pc has been 
detected: such a cloud would have covered most of the 
HFPs in our sample. Thus, our results are possibly indi- 
cating that the filling factor of such clouds is low, as it 
could be the case in 4C 12.50. 

Finally, the great majority of the fast outflows detected 
in both pc- and kpc-scale radio galaxies have an H I opti- 
cal depth of only r < 0.006 (Morganti et al. 2005). Our 
observations are not sensitive enough to reach such low 
values, therefore the lack of detection of broad absorp- 
tion (except for the case of OQ 208 where the continuum 
is indeed strong enough) indicates that H I outflows are 
characterized by low optical depths. 

6. Conclusions 

An Hi absorption search has been carried out with the 
WSRT for a sample of 6 HFP radio galaxies. We con- 
firm the detection of Hi in absorption in 2 galaxies. One 
source, J0111+3906, is characterised by a line width of ~ 
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Fig. 5. Cartoon of a possible orientation of the HFP sources 
and circumnuclear torus. The toroidal structure is perpendic- 
ular to the source axis. The torus model is dependent on the 
opening angle 9 and the viewing angle a. The scale shown is 
approximate: the inner radius of the torus and the sizes of the 
denser clumps are not known. 



100 km s" 1 and a high optical depth of r = 0.44. In the 
other source (OQ 208), the line profile is very broad (~ 
1800 km s _1 ), blue-shifted and shallow, with a maximum 
optical depth of r = 0.005. In the remaining 4 galaxies no 
evidence of Hi absorption has been detected. 

Although this result does not seem to follow the inverse 
correlation found (Pilhstrom et al. 2003) between the col- 
umn density and the linear sizes, it can be explained by 
orientation effects in a torus scenario, in which our line 
of sight intersects the torus in its inner region where the 
low optical depth is due to high spin and kinetic tempera- 
ture. Since these 4 sources have faint flux densities, optical 
depths < 0.01 are not detectable due to sensitivity limita- 
tions. Therefore, the Hi absorption in an object with the 
same spectral characteristics of OQ 208 but with a fainter 
flux density, cannot be detected by our observations. 

Although HFP sources do not seem to follow the corre- 
lation between H I column density and linear size found for 
CSS/GPS sources, this does not imply that we are look- 
ing at a different class of objects, instead of the youngest 
tail of a radio source population. Our results suggest that 
on linear scales smaller than few tens of parsecs, the Hi 
column density is much lower than one would have ex- 
pected on the basis of the work of Pilhstrom et al. (2003). 
As a consequence, compact and rather faint (due to self- 
absorption) sources, such as these HFPs, are not the most 
suitable class of objects to investigate the H I absorption 
on such a small scale. 
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